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ABSTRACT 

High  spatial  resolution  observations  of  the  Ha-emitting  wind  structure  associated  with  the  luminous  blue  variable 
star  P  Cygni  were  obtained  with  the  Navy  Prototype  Optical  Interferometer.  These  observations  represent  the  most 
comprehensive  interferometric  data  set  on  P  Cyg  to  date.  We  demonstrate  how  the  apparent  size  of  the  Ha-emitting 
region  of  the  wind  structure  of  P  Cyg  compares  between  the  2005,  2007,  and  2008  observing  seasons  and  how 
this  relates  to  the  Ha  line  spectroscopy.  Using  the  data  sets  from  2005,  2007,  and  2008  observing  seasons,  we  fit  a 
circularly  symmetric  Gaussian  model  to  the  interferometric  signature  from  the  Ha-emitting  wind  structure  of  P  Cyg. 
Based  on  our  results,  we  conclude  that  the  radial  extent  of  the  Ha-emitting  wind  structure  around  P  Cyg  is  stable  at  the 
10%  level.  We  also  show  how  the  radial  distribution  of  the  Ha  flux  from  the  wind  structure  deviates  from  a  Gaussian 
shape,  whereas  a  two-component  Gaussian  model  is  sufficient  to  fully  describe  the  Ha-emitting  region  around  P  Cyg. 

Key  words:  stars:  individual  (P  Cyg)  -  stars:  winds,  outflows  -  techniques:  interferometric 
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1.  INTRODUCTION 

The  luminous  blue  variable  (LBV)  star  P  Cyg  (HD  193237, 
B2pe)  is  an  unusual  star  with  a  unique  stellar  wind  structure. 
LB  Vs  are  evolved,  very  luminous,  massive  supergiant  stars  that 
show  some  type  of  instability.  As  the  likely  progenitors  of 
Wolf-Rayet  stars  (Crowther  2007)  LB  Vs  may  provide  insight 
into  the  ultimate  end  stage  of  these  stars,  the  core-collapse 
supernovae.  P  Cyg  was  discovered  in  1600  after  a  violent  mass- 
loss  event  caused  it  to  quickly  brighten  to  a  third  magnitude  star 
(de  Groot  1969).  The  violent  mass-loss  events  that  spur  these 
eruptions  are  not  fully  understood,  but  a  number  of  competing 
models  including  single  and  binary  star  scenarios  attempt  to 
explain  these  eruptions  (Humphreys  &  Davidson  1994).  With 
a  mass  loss  between  2  x  10~5  and  4  x  10~4  M0  yr”1  (van 
Blerkom  1978;  Abbott  et  al.  1980;  Leitherer  &  Zickgraf  1987), 
P  Cyg  displays  characteristics  that  are  suggestive  of  energetic 
mass  outflows.  It  is  for  this  reason  that  P  Cyg  has  a  unique  stellar 
wind  structure  with  possible  high-density  regions  in  motion 
within  the  wind  structure.  Despite  numerous  studies  published 
in  the  literature  that  discuss  P  Cyg  or  other  LBVs,  there  is  still 
need  for  observations  that  can  spatially  resolve  the  circumstellar 
region  around  P  Cyg. 

The  spectrum  of  P  Cyg  reveals  emission  line  profiles  with 
their  archetypical  shapes,  which  are  formed  by  the  circumstellar 
matter  that  has  been  shed  from  the  central  star.  These  profiles 
are  commonly  identified  by  an  emission  component  to  the 
red  side  of  an  absorption  line.  The  origin  of  this  Doppler 
shift  can  be  explained  with  a  spherically  symmetric  outflowing 
wind  in  which  the  velocity  increases  with  radial  distance  (see, 
e.g.,  Lamers  &  Cassinelli  1999,  Section  2.2).  The  absorption 
component  corresponds  to  the  continuum  light  absorbed  by  the 
wind  structure  directly  between  the  observer  and  the  central 
star.  The  emission  component  on  the  other  hand  originates  in 
the  spherically  symmetric  halo  around  the  central  star.  Because 
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the  same  characteristics  apply  to  the  Ha  emission  line,  which 
is  one  of  the  strongest  emission  lines  in  the  spectrum  of  P  Cyg, 
this  makes  the  Ha  emission  line  an  excellent  probe  of  the  wind 
region  around  P  Cyg. 

P  Cyg  currently  has  an  estimated  mass  of  30  ±  10  Me 
(Lamers  et  al.  1983,  1985).  At  the  start  of  its  lifetime,  P  Cyg 
possibly  had  an  initial  mass  of  50  M0  (de  Groot  &  Lamers 
1992),  but  has  shed  much  of  its  mass  because  of  its  violent 
history  of  mass-loss  events.  The  effective  temperature  and  the 
stellar  diameter  of  P  Cyg  are  estimated  to  be  Teff  =  19,300  ± 
2000  K  and  76  ±  14  A’0,  respectively  (Lamers  et  al.  1983). 
However,  the  Ha -emitting  region  extends  radially  much  farther 
than  the  central  star  allowing  us  to  resolve  this  region  using  long- 
baseline  interferometric  techniques,  even  though  the  distance  is 
estimated  to  be  in  the  range  from  1 .7  ±  0. 1  kpc  to  1 .8  ±  0. 1  kpc 
(Lamers  et  al.  1983;  Najarro  et  al.  1997). 

The  circumstellar  region  of  P  Cyg  has  been  spatially  resolved 
in  the  past  using  radio  and  optical  interferometry,  as  well  as 
direct  imaging  with  adaptive  optics  (AO).  Radio  interferometric 
observations  detect  the  nebula  around  P  Cyg  at  the  angular 
scales  from  ~50  mas  to  almost  an  arcminute  (Skinner  et  al. 
1997,  1998),  whereas  optical  long-baseline  interferometry  is 
sensitive  to  structures  in  the  nebula  of  P  Cyg  that  are  more  than 
an  order  of  magnitude  smaller.  Therefore,  optical  interferometry 
provides  a  unique  window  of  opportunity  to  resolve  the  inner 
portion  of  the  wind  structure. 

Most  relevant  to  our  study  are  the  results  that  probed  the 
Ha-emitting  region  of  P  Cyg.  For  example,  based  on  near- 
diffraction-limited  observations  obtained  using  the  AO  system 
on  a  1.52  m  telescope,  Chesneau  et  al.  (2000)  not  only  resolved 
the  outer  Ha-emitting  region  of  the  extended  envelope,  but 
detected  signatures  of  clumping.  Although,  the  angular  scales 
sampled  with  a  1.52  m  telescope  were  quite  large,  in  the  region 
of  600  mas,  the  angular  resolution  was  limited  to  50  mas.  On 
much  smaller  angular  scales,  Vakili  et  al.  (1997)  using  a  single 
17.7  m  baseline  of  the  GUT  interferometer  have  resolved  the 
inner  Ha-emitting  envelope  of  P  Cyg  yielding  an  angular  size 
of  5.52  ±  0.47  mas  and  reported  a  marginal  detection  at  the 
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He  i  6678  line.  The  spatial  scales  probed  by  our  study  are 
similar  to  those  sampled  by  Vakili  et  al.  (1997);  however,  the 
interferometric  observations  at  the  Ha  line  obtained  in  our  study 
were  acquired  using  1 1  unique  interferometric  baselines  ranging 
in  length  from  15.8  to  79.3  m,  thus  resulting  in  much  higher 
angular  resolution  measurements. 

2.  OBSERVATIONS 
2.1.  Interferometry 

High  angular  resolution  long-baseline  interferometric  obser¬ 
vations  were  obtained  on  P  Cyg  over  37  nights  in  2005, 2007,  and 
2008  using  the  Navy  Prototype  Optical  Interferometer  (NPOI) 
located  near  Flagstaff,  Arizona.  The  NPOI  is  a  long-baseline 
optical  interferometer  which  consists  of  six  50  cm  siderostat 
mirrors  that  send  12.5  cm  diameter  light  beams  down  the  opti¬ 
cal  feed  system  consisting  of  15  cm  diameter  pipes  to  an  optics 
lab  where  the  interference  (i.e.,  a  fringe)  between  the  beams  is 
recorded  (Armstrong  et  al.  1998).  Observations  covering  a  wide 
spectral  band  can  be  taken  over  multiple  baselines  (currently 
ranging  from  19  to  80  m  in  length)  and  are  geometrically  com¬ 
pensated  and  modulated  using  vacuum  delay  lines  to  guarantee 
that  the  entire  band  is  being  observed  simultaneously. 

The  observing  process  involves  measuring  the  fringe  contrast 
over  15  spectral  channels  every  2  ms.  Each  2  ms  frame  is 
then  used  to  obtain  a  squared  visibility  (V2)  measure  (i.e., 
a  V2  value  measures  the  fringe  contrast).  Afterward,  these 
values  are  incoherently  averaged  over  1  s  intervals  to  create 
1  s  data  points  (Hummel  et  al.  2003).  The  final  observations 
consist  of  30  s  of  integration  (also  known  as  a  “scan”)  during 
which  the  squared  visibilities  are  measured  simultaneously  at 
three  baselines  per  output  from  the  beam  combiner  (with  two 
outputs  being  used  simultaneously)  and  15  spectral  channels 
covering  the  wavelengths  between  560  and  870  nm.  The  typical 
observational  sequence  consists  of  a  pair  of  coherent  and 
incoherent  (off  the  fringe)  scans  on  a  target  followed  by  a 
pair  of  scans  on  a  calibrator.  The  incoherent  scans  are  used  to 
estimate  the  additive  bias  terms  affecting  the  squared  visibility 
measures.  The  process  is  repeated  during  a  night  for  as  long  as 
the  fringe  packet  can  be  detected  on  the  target  star.  Typically, 
the  observations  stop  when  the  optical  path  difference  between 
the  different  siderostats  cannot  be  compensated  anymore. 

The  Ha -emitting  sources  form  a  specific  class  of  targets  for 
NPOI  since  the  entire  emission  line  ends  up  only  in  one  out 
of  the  15  spectral  channels.  This  characteristic  allows  one  to 
calibrate  the  V2  values  from  the  channel  containing  the  Ha 
emission  using  the  V2  values  from  the  neighboring  continuum 
channels  (Tycner  et  al.  2003).  Although  during  this  calibration 
process  one  needs  to  assume  that  the  diameter  of  the  central  star 
at  the  continuum  wavelengths  is  known,  the  final  results  based 
on  the  Ha  emission  are  only  weakly  dependent  on  the  adopted 
value  for  the  angular  diameter  of  the  central  star  (we  explore  this 
in  more  detail  in  Section  3.2).  For  the  purpose  of  the  calibration, 
we  adopted  a  uniform  disk  (UD)  diameter  of  0.2  mas  to  represent 
the  central  star  at  the  continuum  wavelengths.  We  obtain  this 
diameter  based  on  a  distance  estimate  of  1.8  ±  0.1  kpc  and  a 
central  star  radius  of  76  ±  14  Re  (Lamers  et  al.  1983). 

Applying  the  procedures  described  above  to  the  entire  obser¬ 
vational  data  set  we  acquired  for  P  Cyg,  yielded  a  total  of  1534 
Ha  squared  visibilities  on  1 1  unique  baselines  (out  of  a  max¬ 
imum  of  15  unique  baselines  for  a  six-station  configuration). 
Table  1  lists  the  number  of  unique  Ha  V2  measures  we  obtained 
for  each  of  the  nights  in  our  observing  campaign.  Although  the 


Table  1 

NPOI  Observing  log  of  P  Cyg 


UT  Date 

No.  of  Ha  squared 

Visibilities 

2005 Jun  2 

48 

2005  Jun  3 

19 

2005  Jun  4 

48 

2005  Jun  6 

18 

2005  Jun  8 

20 

2005  Jun  9 

84 

2005  Jun  1 1 

19 

2005 Jun  13 

50 

2005  Jun  16 

90 

2007  Jul  1 

20 

2007  Jul  2 

20 

2007  Jul  4 

20 

2007  Jul  12 

15 

2007  Aug  8 

24 

2007  Aug  9 

16 

2007  Aug  10 

66 

2007  Aug  1 1 

56 

2007  Aug  12 

28 

2007  Aug  18 

72 

2007  Aug  19 

31 

2007  Aug  20 

82 

2007  Aug  21 

33 

2007  Aug  22 

92 

2007  Aug  23 

17 

2007  Aug  24 

12 

2007  Aug  25 

41 

2007  Aug  27 

12 

2007  Aug  29 

12 

2007  Sep  5 

36 

2008  Jun  27 

16 

2008  Jun  28 

40 

2008  Jun  30 

36 

2008  Jul  2 

28 

2008  Jul  3 

60 

2008  Jul  7 

64 

2008  Jul  8 

89 

2008  Jul  9 

100 

Total 

1534 

observational  setup  for  each  of  the  three  observing  seasons  was 
slightly  different,  and  therefore  a  range  of  different  baselines 
was  sampled  each  year,  in  total  we  have  obtained  observations 
at  baselines  ranging  from  the  shortest  (15.8  m)  to  the  longest 
(79.3  m)  possible  at  the  time  of  our  observations.  More  specif¬ 
ically,  the  2005  and  2008  seasons  contained  baselines  between 
18.9  and  53.2  m,  and  the  2007  season  contained  baselines  from 
15.8  to  79.3  m  in  length.  Figure  1  shows  the  sampling  of  the 
(m,  i;)-plane  in  the  Ha  channel  for  P  Cyg  based  on  all  the  avail¬ 
able  nights  in  2005,  2007,  and  2008.  The  “arcs”  seen  in  the  plot 
are  a  result  of  the  changing  projection  of  the  baseline  vector  on 
the  sky  during  the  night  due  to  Earth’s  rotation.  Three  seasons 
of  squared  visibility  data  sampled  across  the  (n,  i>)-plane  are 
shown  in  Figure  2  plotted  as  a  function  of  radial  spatial  fre¬ 
quency  (V m2  +  v2).  The  calibrated  Ha  V2  measures  are  listed 
in  Table  2. 

In  addition  to  the  V2  measures,  we  have  also  obtained  clo¬ 
sure  phase  quantities  (sum  of  complex  phases  on  three  base¬ 
lines  forming  a  triangle)  on  four  unique  baseline  triangles.  By 
fitting  a  quadratic  function  to  the  closure  phases  from  the  con¬ 
tinuum  channels  (at  each  scan  and  closure  phase  configura¬ 
tion  separately)  and  then  removing  this  quadratic  trend  from  all 
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Table  2 

Calibrated  Ha  Squared  Visibilities  of  P  Cyg 


Julian  Date 
(JD  -  2,450,000) 

Spatial  Frequency  u 
(106  cycles  rad-1) 

Spatial  Frequency  v 
(106  cycles  rad-1) 

vL 

Baseline 

3523.863 

-18.364 

19.654 

0.666  ±  0.028 

AC-AE 

3523.863 

-11.184 

-32.672 

0.623  ±  0.040 

AC-AN 

3523.863 

-7.180 

52.326 

0.615  ±  0.017 

AE-AN 

3523.863 

-37.854 

29.298 

0.567  ±  0.023 

AC-E6 

3523.863 

67.722 

-25.810 

0.512  ±0.026 

AW-E6 

3523.863 

29.868 

3.488 

0.671  ±  0.045 

AC-AW 

3523.891 

-21.636 

17.482 

0.705  ±  0.023 

AC-AE 

3523.891 

-8.106 

-33.720 

0.633  ±  0.039 

AC-AN 

3523.891 

-13.530 

51.202 

0.515  ±  0.023 

AE-AN 

3523.891 

-43.328 

24.890 

0.526  ±  0.025 

AC-E6 

3523.891 

74.493 

-18.088 

0.501  ±  0.033 

AW-E6 

(This  table  is  available  in  its  entirety  in  machine-readable  and  Virtual  Observatory  (VO)  forms  in  the  online  journal.  A 
portion  is  shown  here  for  guidance  regarding  its  form  and  content.) 


Figure  1.  Sampling  of  the  (m,  u)-plane  in  the  Ha  channel  for  P  Cyg  on  multiple 
baselines  ranging  in  length  from  15.8  to  79.3  m  obtained  on  37  nights  of 
observation  in  2005,  2007,  and  2008  (open  circles).  Possible  coverage  at  each 
baseline  is  also  shown  from  6  hr  east  of  meridian  (dotted  line)  to  6  hr  west  of 
meridian  (solid  line). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

channels,  including  the  Ha  channel,  we  were  able  to  obtain  the 
differential  closure  phases  for  the  Ha-emitting  region.  Because 
we  expect  the  central  star  to  be  mostly  unresolved  at  the  con¬ 
tinuum  channels  and  if  we  assume  that  the  continuum  source  is 
point  symmetric,  then  we  expect  the  closure  phases  to  cluster 
around  a  zero  value.  This  allows  us  to  inspect  the  Ha  closure 
phases  with  respect  to  the  mean  zero  closure  phase.  In  order  to 
reduce  the  random  noise  present  in  data  associated  with  a  single 
scan,  we  calculated  the  weighted  mean  closure  phases  for  each 
season  and  each  unique  configuration.  In  the  2007  season,  one 
of  the  baselines  (AE-AN)  was  sampled  at  two  independent  out¬ 
puts  from  the  beam  combiner,  and  this  allowed  us  to  calculate 
two  closure  phase  quantities  containing  this  baseline  either  from 
one  output  beam  or  the  other.  Because  the  observations  for  the 
other  baselines  were  identical,  we  decided  not  to  average  these 
quantities,  and  instead  treat  them  separately.  Figure  3  shows  the 


Figure  2.  Calibrated  squared  visibilities  from  the  Ha  channel  of  P  Cyg  obtained 
on  1 1  unique  baselines  from  37  nights  of  observation  in  2005, 2007,  and  2008.  A 
signature  of  a  central  star  represented  by  a  UD  model  with  an  angular  diameter 
of  0.2  mas  is  shown  (dashed  line)  along  with  three  models  representing  the 
Ha -emitting  envelope.  The  envelope  component  is  modeled  with  a  UD  (dotted 
line),  a  single-component  Gaussian  (dash-dotted  line),  and  a  two-component 
Gaussian  model  (solid  line). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

weighted  mean  closure  phases  for  all  the  configurations  with 
the  quadratic  trend  based  on  the  continuum  channels  already 
removed. 

2.2.  Spectroscopy 

The  Ha  emission  line,  as  detected  by  the  NPOI,  falls  on  a 
spectral  channel  that  has  a  spectral  width  of  15  nm.  Therefore, 
the  interferometric  data  obtained  on  P  Cyg  do  not  provide 
sufficient  spectroscopic  information  to  allow  us  to  determine  the 
strength  or  the  shape  of  the  Ha  emission  line.  For  this  reason,  we 
obtained  complementary  spectroscopy  of  P  Cyg  using  Lowell 
Observatory’s  John  S.  Hall  1.1m  telescope  located  at  the  same 
observing  site  as  the  NPOI.  The  telescope  was  outfitted  with 
the  solar-stellar  spectrograph  (SSS),  which  is  a  fiber-fed  echelle 
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Figure  3.  Weighted  mean  closure  phases  (with  a  quadratic  trend  removed  from  the  continuum  channels)  obtained  at  each  possible  baseline  triangle  obtained  in  2005, 
2007,  and  2008  observing  seasons.  The  signal  from  the  channel  containing  the  Ha  emission  is  marked  with  red  squares.  The  closure  phases  containing  the  AC-AE-AN 
and  E6-AE-AN  stations  in  the  2007  season  are  plotted  in  two  panels  since  one  of  the  baselines  (AE-AN)  was  observed  at  two  independent  output  beams  from  the 
beam  combiner  (indicated  in  the  square  brackets  with  an  SI  or  S2  designator),  thereby  allowing  calculation  of  two  closure  phase  quantities.  The  error  bars  are  based 
on  the  uncertainty  of  the  weighted  average. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


spectrograph  (Hall  &  Lockwood  1995).  Spectroscopic  data  on 
P  Cyg  were  obtained  during  the  period  from  2005  April  to  2007 
December,  with  the  individual  nights  listed  in  Table  3.  The  SSS 
instrument  produces  spectra  in  the  Ha  region  with  a  resolving 
power  of  ~10,000.  A  typical  Ha  line  profile  normalized  with 
respect  to  the  continuum  level  of  P  Cyg  is  shown  in  Figure  4.  The 
equivalent  width  (EW)  measurements  based  on  the  individual 
spectra  are  also  plotted  in  the  upper  panel  of  Figure  5. 


3.  ANALYSIS 

3.1.  Interferometric  Results 

There  are  two  signatures  in  the  squared  visibilities  from  the 
spectral  channel  containing  the  Ha  emission  line  shown  in 
Figure  2.  One  is  from  the  central  star  and  the  other  is  from 
the  Ha-emitting  circumstellar  region.  To  model  these  data,  we 
use  the  same  two-component  model  as  was  used  by  Tycner  et  al. 
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Table  3 

Spectroscopic  Variability  for  P  Cyg 


UT  Date 

-EW  (A) 

%  Change3 

^Tnax 

Fc 

%  Change3 

2005  Apr  1 

82.1 

3.5 

20.8 

6.7 

2005  May  19 

74.9 

5.7 

18.7 

4.0 

2005  Jun  28 

81.0 

2.0 

20.8 

6.7 

2005  Jul  16 

80.3 

1.1 

20.8 

6.9 

2005  Aug  22 

79.7 

0.5 

19.7 

1.1 

2005  Sep  16 

78.8 

0.8 

19.1 

2.0 

2005  Sep  17 

81.3 

2.4 

19.6 

0.6 

2005  Sep  27 

81.0 

2.1 

19.5 

0.1 

2005  Oct  1 1 

81.0 

2.1 

18.9 

3.1 

2005  Nov  8 

73.6 

7.3 

17.0 

12.9 

2005  Mean 

79.4  ±  0.9 

19.5  ±  0.4 

2006  Jun  10 

81.4 

7.3 

19.7 

7.3 

2006  Aug  3 

70.1 

7.7 

16.6 

9.4 

2006  Sep  6 

75.3 

0.8 

17.9 

2.7 

2006  Oct  3 

70.3 

7.3 

17.2 

6.2 

2006  Oct  28 

78.1 

2.9 

19.0 

3.3 

2006  Oct  31 

77.4 

2.0 

19.0 

3.4 

2006  Nov  1 

78.6 

3.5 

19.2 

4.4 

2006  Mean 

75.9  ±  1.6 

18.4  ±0.4 

2007  May  3 

70.5 

10.8 

17.3 

17.7 

2007  Jun  4 

75.8 

4.1 

18.7 

11.1 

2007  Jun  22 

80.0 

1.2 

19.7 

6.4 

2007  Oct  23 

83.4 

5.5 

22.5 

6.8 

2007  Oct  24 

78.0 

1.3 

21.1 

0.3 

2007  Oct  27 

76.5 

3.2 

20.9 

0.6 

2007  Oct  28 

74.7 

5.6 

20.2 

3.9 

2007  Dec  18 

80.9 

2.4 

22.6 

7.5 

2007  Dec  19 

85.5 

8.2 

23.7 

12.7 

2007  Dec  20 

85.2 

7.8 

23.7 

12.5 

2007  Mean 

79.0  ±  1.6 

21.0  ±0.7 

Note.  3  Changes  calculated  with  respect  to  yearly  mean  EW  and  peak  strength. 


(2006)  of  the  form 

^model  =  [cp^p  +  0  —  Cp)Venv]  ,  (1) 

where  Vp  is  the  visibility  function  representing  the  photosphere 
of  the  central  star,  Venv  is  the  visibility  function  representing 
the  circumstellar  envelope,  and  cp  is  the  fractional  contribution 
from  the  stellar  photosphere  to  the  total  flux  in  the  Ha  channel. 

Based  on  the  results  from  differential  closure  phases  (recall 
Section  2. 1),  which  showed  very  weak  (if  any)  variations  across 
the  Ha  channel,  we  expect  the  Ha-emitting  structure  to  be 
sufficiently  well  described  (at  the  level  of  angular  resolution 
provided  by  our  range  of  baselines)  by  a  model  that  is  both 
point  symmetric  and  concentric  with  the  central  star.  This  means 
that  Vp  and  Venv  in  Equation  (1)  can  be  represented  by  real 
functions  (i.e.,  the  Fourier  transform  of  real  and  even  function 
results  in  a  real  function).  To  model  the  stellar  photosphere 
we  use  a  circularly  symmetric  UD  brightness  distribution.  The 
normalized  visibility  amplitude  for  a  UD  can  be  written  as 


Vp(u,  v)  =  2 


Ji(jtas ) 
Jtas 


(2) 


where  J  \  denotes  the  Bessel  function  of  the  first  kind  and  first 
order,  a  is  the  angular  diameter,  and  s  is  the  radial  spatial 
frequency  (i.e.,  s  =  V u1  +  v2).  Because  the  interferometric 
observations  of  P  Cyg  were  obtained  at  sufficiently  high  spatial 
frequencies  so  that  the  Ha-emitting  region  was  fully  resolved, 
we  can  compare  different  models  of  brightness  distribution  for 


Wavelength  (nm) 


Figure  4.  Representative  Ha  line  profile  of  P  Cyg  obtained  on  2007  June  22. 
The  EW  of  the  emission  line  is  —80.0  A. 


the  envelope  component.  We  choose  to  model  our  data  with  a 
UD  and  a  Gaussian  distribution  (GD),  where  in  the  former  Venv 
has  the  same  functional  form  as  Equation  (2)  except  that  a  is 
replaced  with  Pud,  and  in  the  latter  case 


yGD 

env 


=  exp 


(n9GDs)2 

4  In  2 


(3) 


where  Pgd  corresponds  to  the  FWHM  of  the  Gaussian  model. 

The  envelope  component  (Venv)  in  Equation  (1),  represented 
by  either  a  UD  or  a  Gaussian  radial  distribution,  was  assumed 
to  be  circularly  symmetric.  This  choice  was  based  on  a  test 
performed  on  the  observational  data  set  being  divided  into  18 
subsets  based  on  the  positional  angles  (P.A.s)  covered  by  the 
observations.  Figure  6  shows  the  variation  of  the  fitted  diameter 
of  a  Gaussian  model  as  a  function  of  PA.  Although  a  hint  of 
variation  is  present  in  the  data,  more  than  68%  (based  on  ±ler 
expectation)  of  points  fall  within  5.5%  of  the  mean  value,  and 
thus  we  conclude  that  some  of  that  variation  might  be  solely 
due  to  the  fact  that  the  observations  were  acquired  over  many 
observing  seasons  and  intrinsic  variability  at  the  5%  level  cannot 
be  ruled  out.  A  similar  test  performed  on  much  smaller  data  set 
from  only  one  season  does  not  reveal  any  variation  beyond  those 
expected  based  on  the  uncertainties  associated  with  the  fitted 
model  parameters.  Therefore,  for  the  purpose  of  describing  the 
geometrical  structure  of  the  Ha-emitting  region  of  P  Cyg,  we 
decided  to  concentrate  only  on  circularly  symmetric  models 
fitted  to  the  V 2  data  from  the  individual  observing  seasons,  as 
well  as  to  the  entire  data  set. 

Table  4  summarizes  both  UD  and  GD  model  results  along 
with  their  corresponding  reduced  /2  values  (x2),  which  were 
used  to  assess  the  goodness  of  fit  of  each  model.  Based  on  the  x2 
values  obtained  from  model  fits  to  the  entire  data  set,  it  is  clearly 
evident  that  the  Gaussian  model  (with  x2  =  1-9)  produced  a 
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Julian  Day  (2440000+) 


Julian  Day  (2440000+) 

Figure  5.  Upper  panel:  P  Cyg  H»  EWs  from  1988  July  to  2007  December  with  data  from  Markova  et  al.  (2001a)  (triangles)  and  this  study  (squares)  shown.  The 
vertical  lines  mark  the  periods  when  interferometric  observations  were  obtained  at  the  NPOI  (dashed  lines),  and  the  single  night  of  interferometric  observation  of  P 
Cyg  reported  by  Vakili  et  al.  (1997)  (dotted  line).  Lower  panel:  the  AAVSO  Johnson  V-band  photoelectric  light  curve  of  P  Cyg  for  the  period  2005-2009,  obtained 
based  on  measurements  relative  to  the  comparison  star  HD  188892  ( V  =  4.936,  B  —  V  =  —0.086).  Only  data  where  the  measured  check  star  magnitude  (HD  193369; 
V  =  5.573)  lies  within  ±0.04  mag  of  the  known  value  are  plotted.  Typical  photometric  errors  per  point  are  5—10  mmag.  The  smoothed  moving  average  (with  a  1000 
d  window)  is  also  shown  (solid  line)  to  guide  the  eye  for  possible  long-term  variability. 


Table  4 

Best-fit  Model  Parameters  for  P  Cyg 


Model 

Season 

Diameter  (mas) 

Second  Diameter  (mas) 

k\  or  l\ 

Cp 

Xv 

Uniform  disk3 

2005 

9.41  ±  2.50 

0.84  ±  0.03 

Uniform  disk 

2005 

6.74  ±  0.05 

0.78  ±  0.01 

2.6 

Gaussian 

2005 

3.97  ±  0.05 

0.75  ±  0.01 

1.5 

Double-Gaussian 

2005 

8.7  ±  5.3 

2.08  ±  0.15 

0.41  ±  0.03 

0.70  ±  0.01 

1.2 

Uniform  disk3 

2007 

10.17  ±  0.32 

0.84  ±  0.01 

Uniform  disk 

2007 

7.37  ±  0.03 

0.79  ±  0.01 

3.7 

Gaussian 

2007 

4.37  ±  0.03 

0.76  ±  0.01 

2.4 

Double-Gaussian 

2007 

5.33  ±  0.15 

1.53  ±  0.15 

0.60  ±  0.02 

0.70  ±  0.01 

1.7 

Uniform  disk3 

2008 

8.40  ±  1.50 

0.79  ±  0.04 

Uniform  disk 

2008 

7.06  ±  0.06 

0.78  ±  0.01 

1.5 

Gaussian 

2008 

4.37  ±  0.07 

0.75  ±  0.01 

1.3 

Double-Gaussianb 

2008 

10  ±  9 

2.17  ±  0.07 

0.50  ±  0.02 

0.70 

1.1 

Uniform  disk3 

All 

9.83  ±  0.26 

0.84  ±  0.01 

Uniform  disk 

All 

7.21  ±  0.03 

0.79  ±  0.01 

2.8 

Gaussian 

All 

4.29  ±  0.03 

0.76  ±  0.01 

1.9 

Double-Gaussian 

All 

5.64  ±  0.17 

1.80  ±  0.11 

0.57  ±  0.02 

0.70  ±  0.01 

1.5 

Uniform  disk  +  Gaussian 

All 

3.06  ±  0.15 

5.46  ±  0.16 

0.36  ±  0.02 

0.72  ±  0.01 

1.5 

Notes. 

Ll  Models  fitted  to  data  for  baselines  up  to  18.9  m. 

b  There  were  insufficient  data  at  the  longest  baselines  to  fit  for  four  model  parameters  and  instead  a  fit  was  obtained  with  cp  fixed  at  0.70. 


significantly  better  fit  to  the  data  than  a  circularly  symmetric 
UD  model  (with  xl  —  2.8).  However,  neither  the  UD  nor  the 
GD  model  fully  reproduces  the  trend  seen  in  the  data  obtained 
at  high  spatial  frequencies  as  can  be  seen  in  Figure  2.  Although, 


the  Gaussian  model  does  produce  a  better  fit,  to  fully  account 
for  the  observational  signature  at  the  high  spatial  frequencies  we 
require  a  model  with  an  extra  degree(s)  of  freedom.  A  natural 
extension  to  a  Gaussian  model  is  a  two-component  Gaussian 
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PA  (deg) 


Figure  6.  Diameter  of  a  Gaussian  model  fitted  to  subsets  of  squared  visibility 
data  divided  based  on  ranges  of  P.A.  on  the  sky.  Each  point  represents  a  10° 
range  in  P.A.  The  mean  diameter  of  4.29  d=  0.03  mas  based  on  a  model  fit  to  all 
data  (dashed  line)  and  5.5%  variation  around  that  value  (dotted  lines)  are  also 
shown. 


model  of  the  form 


y2GD 


=  k\  exp 


(JT0\ S)2 

4  In  2 


+  (1  —  £i)  exp 


(jtd2s)2 

4  In  2 


(4) 


where  0\  and  02  correspond  to  the  FWHM  values  of  the  two 
Gaussian  components,  and  k\  is  the  fractional  contribution 
from  the  first  Gaussian.  Similar  two-component  Gaussian  fits 
have  been  applied  to  other  stars  by  Hofmann  et  al.  (2002) 
and  Preibisch  et  al.  (2003).  Combining  Equations  (1),  (2), 
and  (4),  and  fitting  the  two-component  Gaussian  model  to  the 
observations  from  all  three  observing  seasons,  results  in  angular 
diameters  of  0\  =  5.64  ±  0.17  mas  and  02  —  1.80  ±  0.11  mas, 
with  a  fractional  contribution  from  the  stellar  photosphere, 
cp  =  0.70  ±  0.01  and  k\  =  0.57  ±  0.02  (i.e.,  with  the  two 
Gaussian  components  contributing  approximately  equally  to  the 
net  Ha  emission).  The  best-fit  two-component  Gaussian  model 
had  the  lowest  reduced  x  2  value  (yielding  /  “  of  1 .5  based  on  the 
entire  data  set)  out  of  all  the  models  considered  in  our  analysis. 
The  model  curve  corresponding  to  the  two-component  Gaussian 
model  is  also  shown  in  Figure  2. 

The  success  of  the  two-component  Gaussian  model  to  repro¬ 
duce  the  observational  signature  from  the  Ha  channel  is  mostly 
due  to  the  higher  degree  of  freedom  along  the  radial  direction. 
In  other  words,  one  can,  in  principle,  approximate  any  mono- 
tonically  decreasing  function  as  a  sum  of  a  large  number  of 
Gaussian  functions.  In  fact,  any  other  monotonically  decreasing 
function  can  also  be  used  instead  of  a  Gaussian.  For  example,  a 
sum  of  a  GD  and  a  UD  component  of  the  form 


yUD+GD 
'  env 


7i(jt0ud^) 

l\Z - 

7TO\jdS 


+  (1  —  / 1)  exp 


(tt0gds)2 

4  In  2 


where  6>ud  and  @gd  correspond  to  the  angular  diameters  of 
the  UD  and  Gaussian  components,  respectively,  and  /]  is  the 


fractional  contribution  from  the  UD  and  yields  equally  good  fit 
to  the  data  as  the  two-component  Gaussian  model  (resulting  in 
a  reduced  x2  of  1.5).  The  best-fit  angular  diameters  are  in  this 
case  0\jd  =  3.06  ±  0.15  mas  and  6>gd  =  5.46  ±  0.16  mas,  with  a 
fractional  contribution  from  the  stellar  photosphere,  cp,  of  0.72 
and  /[  =  0.36. 


3.2.  Central  Star  Diameter 

A  central  star  diameter  of  0.2  mas  was  assumed  for  our  model 
(recall  Section  2.1).  Because  varying  the  adopted  diameter 
for  the  central  star  affects  the  calculated  size  of  the  Ha- 
emitting  region,  we  test  the  dependency  of  the  angular  diameter 
of  the  Ha-emitting  wind  structure  surrounding  P  Cyg  on 
the  angular  diameter  of  the  photospheric  component.  This 
is  done  by  fitting  Equation  (1)  to  the  observations  with  the 
envelope  component  represented  by  a  Gaussian  model  and 
the  photospheric  component  with  a  UD  model  for  a  range  of 
adopted  stellar  diameters.  Using  the  observations  from  the  2007 
observing  season  we  find  that  the  assumed  central  star  diameter 
of  P  Cyg  at  most  has  only  a  1%  effect  on  the  derived  angular 
size  of  the  wind  structure  when  the  adopted  angular  diameter  of 
the  central  star  is  varied  from  0  to  1  mas.  Therefore,  an  angular 
size  of  the  continuum-emitting  region  that  is  underestimated  by 
a  factor  of  5  results  only  in  a  1%  effect  on  our  best-fit  angular 
diameter  of  the  Ha-emitting  region. 

3.3.  Spectroscopic  Results 

To  test  for  any  possible  link  between  the  changes  in  the 
strength  of  the  Ha  emission  line  and  any  variations  in  the  size  of 
the  Ha-emitting  region,  we  examined  the  spectroscopic  data  for 
any  variability  on  a  timescale  comparable  to  our  interferometric 
observations.  We  found  that  there  was  EW  variability  with  an 
amplitude  <10  A  over  a  timescale  of  ~2  months,  which  is 
consistent  with  the  short-term  variability  seen  by  Markova  et  al. 
(2001b).  However,  the  yearly  mean  Ha  EW  and  the  yearly  mean 
Ha  peak  signal  showed  no  variation  from  year  to  year  during  our 
observing  period  suggesting  relative  stability  in  EW  measure  on 
a  year-to-year  timescale  (see  Table  3). 

While  our  mean  spectroscopic  data  suggest  stability  on  a 
year-to-year  scale,  we  cannot  rule  out  variability  of  Ha  emission 
if  such  changes  are  correlated  with  changes  in  the  continuum 
level  as  suggested  by  Markova  et  al.  (2001a).  If  the  continuum 
brightened  and  the  Ha  has  increased  in  strength  by  the  same 
fractional  amount,  it  would  appear  stable.  Regardless  of  this 
apparent  mean  stability,  we  detect  high-frequency  variations 
on  timescales  much  shorter  than  one  year.  The  upper  panel  of 
Figure  5  shows  the  Ha  EWs  derived  from  our  spectroscopic 
observations  along  with  the  data  published  by  Markova  et  al. 
(2001a). 

Markova  et  al.  (2001b)  obtained  both  spectroscopic  and 
photometric  observations  of  P  Cyg  and  found  that  the  EW 
measurements  displayed  a  slow  variability,  with  an  amplitude 
of  about  30  A  and  a  duration  of  about  600  days,  and  a  faster 
variability  with  an  amplitude  of  up  to  10  A  and  a  duration  of 
40-60  days.  While  we  were  able  to  observe  a  variability  of  up 
to  10  A  in  amplitude  on  a  short  (month-to-month)  timescale, 
we  did  not  observe  the  long-term  30  A  amplitude  variability. 
If  the  long-term  variability  discussed  in  Markova  et  al.  (2001b) 
was  present  at  the  time  of  their  observations,  it  is  no  longer 
present  in  our  data.  Therefore,  based  on  our  spectroscopic  and 
interferometric  observations,  we  conclude  that  P  Cyg  is  in  a 
quiescent  phase  of  its  lifetime,  and  it  is  not  going  through  any 
drastic  changes. 
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Table  5 

Single  Night  Spectroscopic  Variability  for  P  Cyg 


UT  Date 

-EW  (A) 

%  Change2 

^Tnax 

Fr 

%  Change3 

2007  Dec  20 

85.2 

0.1 

23.7 

0.0 

2007  Dec  20 

85.9 

0.9 

23.8 

0.7 

2007  Dec  20 

84.4 

0.9 

23.5 

0.7 

2007  Dec  20 

85.6 

0.5 

23.8 

0.4 

2007  Dec  20 

85.2 

0.1 

23.7 

0.0 

2007  Dec  20 

85.8 

0.7 

23.8 

0.6 

2007  Dec  20 

86.0 

1.0 

23.8 

0.6 

2007  Dec  20 

83.4 

2.1 

23.4 

1.2 

2007  Dec  20 

85.0 

0.2 

23.6 

0.2 

2007  Dec  20 

85.2 

0.0 

23.7 

0.0 

2007  Dec  20 

83.4 

2.1 

23.4 

1.2 

2007  Dec  20 

84.7 

0.6 

23.6 

0.4 

2007  Dec  20 

85.3 

0.1 

23.6 

0.1 

2007  Dec  20 

85.5 

0.4 

23.7 

0.3 

2007  Dec  20 

85.4 

0.2 

23.7 

0.1 

2007  Dec  20 

85.7 

0.6 

23.8 

0.5 

2007  Dec  20 

86.1 

1.1 

23.8 

0.6 

2007  Dec  20 

85.5 

0.4 

23.7 

0.3 

2007  Dec  20 

85.1 

0.1 

23.6 

0.3 

2007  Dec  20  Mean 

85.2  ±  0.2 

23.67  ±  0.03 

Note.  a  Changes  calculated  with  respect  to  nightly  mean.  All  spectra  were 
obtained  with  180  s  of  integration. 


The  variability  of  the  Ha  EW  in  our  data  set  is  similar 
to  the  10  A  variation  seen  in  the  literature  on  a  timescale 
of  ~2  months.  Based  on  our  analysis  of  the  EW  measures 
listed  in  Table  3,  we  conclude  that  the  rms  variations  are 
less  than  10%  and  typically  closer  to  the  2%-5%  range. 
Therefore,  for  the  purpose  of  our  study  we  assume  that  the 
overall  emission  in  the  Ha  line  throughout  our  interferometric 
run  was  stable  at  the  10%  level.  We  also  tested  for  variability 
in  the  Ha  strength  during  a  single  night  and  found  it  to  be 
virtually  non-variable  as  well  (see  Table  5).  The  variations 
were  within  the  expected  observational  uncertainties  of  ~3%, 
mostly  due  to  uncertainties  associated  with  continuum  level 
normalization. 

4.  DISCUSSION 
4.1.  The  Radial  Intensity  Falloff 

The  deviation  from  a  Gaussian  shape  of  the  spatially  resolved 
wind  structure  of  P  Cyg  is  one  of  our  most  intriguing  results. 
Examining  the  squared  visibilities  in  Figure  2,  it  is  apparent  that 
the  observations  at  high  spatial  frequencies  cannot  be  repre¬ 
sented  by  a  single-component  Gaussian  model.  It  is  worthwhile 
to  mention  that  the  Gaussian  intensity  distribution  is  still  a  closer 
approximation  to  the  radial  intensity  distribution  from  the  wind 
structure  than  the  UD  model  (shown  with  the  dotted  line  in 
Figure  2).  To  our  knowledge,  this  is  the  first  time  the  Ha- 
emitting  region  of  P  Cyg  has  been  observed  to  clearly  deviate 
from  a  Gaussian  shape. 

To  reproduce  our  data  at  the  highest  spatial  frequencies, 
we  use  the  two-component  Gaussian  model  as  a  way  of 
parameterizing  the  spatially  resolved  signature  of  the  Ha- 
emitting  region  of  P  Cyg.  This  suggests  that  the  wind  structure 
is  indeed  more  complex  in  P  Cyg  as  compared  to  Be  stars  (see, 
for  example,  Tycner  et  al.  2008).  Because  P  Cyg  possesses 
radiatively  driven  stellar  wind,  the  velocity  structure  can  be 
described  by  the  /1-law  of  Castor  &  Lamers  (1979),  which 


implies  that  the  wind  is  accelerated  to  its  terminal  velocity  close 
to  the  star  and  then  reaches  a  constant  terminal  velocity.  This  in 
turn  means  that  the  density  structure  of  P  Cyg  will  change  at  two 
different  rates,  one  when  the  wind  is  still  accelerating  and  one 
when  the  wind  has  already  reached  its  terminal  velocity  at  which 
point  the  density  would  be  expected  to  drop  with  radial  distance 
as  r~ 1 .  Although  this  cannot  be  used  as  a  direct  explanation  for 
the  two-component  structure  because  the  Ha  emission  in  P  Cyg 
is  optically  thick  (Leitherer  1988),  by  combining  the  effects 
of  the  density  distribution  with  the  radiative  transfer  effects 
it  should  be  possible  to  use  the  data  presented  in  this  study  to 
directly  constrain  a  numerical  wind  model  based  on  a  /J  velocity 
law.  Lastly,  if  the  two-component  wind  structure  is  dominated 
by  optically  thick  and  thin  effects,  one  could  argue  that  the 
inner  (optically  thick)  region  might  be  better  represented  by  a 
UD  model  and  the  outer  region  by  a  Gaussian  model.  This  would 
be  in  agreement  with  the  results  presented  in  Section  3.1  that 
showed  that  the  observations  can  be  equally  well  represented  by 
a  two-component  UD  and  a  Gaussian  model. 

P  Cyg  has  been  spatially  resolved  in  the  past  using  interfer¬ 
ometric  techniques.  Vakili  et  al.  (1997),  using  a  single  17.7  m 
interferometric  baseline,  resolved  the  extended  envelope  of  P 
Cyg  at  the  Ha  emission  line.  Although  their  data  set  was  limited 
and  did  not  allow  them  to  fit  a  model  more  complex  than  a  UD 
model,  the  angular  diameter  of  the  UD  model  they  obtained  was 
0[  \)  —  5.52  ±  0.47  mas.  They  also  obtained  a  complementary 
Ha  spectrum  on  P  Cyg  on  the  same  night  they  obtained  interfero¬ 
metric  observations  with  a  Ha  peak  signal  of  Fmax/Fcont  =  14.6. 
This  is  somewhat  weaker  than  the  Fmax/Fcont  ~  20  values  we 
obtain  based  on  our  observations  from  2005  to  2007  (recall 
Figure  4),  although  the  effect  of  lower  spectral  resolution  in  the 
spectra  of  Vakili  et  al.  (1997)  cannot  be  ruled  out,  which  could 
explain  the  lower  peak-to-continuum  ratio  seen  in  the  observa¬ 
tion  reported  by  Vakili  et  al.  (1997). 

Although  a  UD  model  is  completely  inconsistent  with  our 
observations  (recall  Figure  2),  in  order  to  compare  our  results 
to  Vakili  et  al.  (1997)  at  effectively  the  same  limiting  spatial 
resolution,  we  fitted  a  UD  model  to  data  from  only  shortest 
baselines  (up  to  18.9  m).  The  best-fit  UD  angular  diameters 
ranged  from  8.4  to  10.2  mas  depending  on  which  season  was 
used  in  the  fit  (see  Table  4).  Therefore,  our  UD  diameter  for  the 
Ha-emitting  wind  structure  is  not  consistent  with  the  reported 
diameter  of  5.52  mas  by  Vakili  et  al.  (1997).  This  discrepancy 
might  be  due  to  the  possibility  that  Vakili  et  al.  (1997)  has 
modeled  the  observational  signal  that  contained  both  the  mostly 
unresolved  central  star  and  the  Ha-emitting  envelope  with  only 
one  UD  diameter.  This  would  be  equivalent  to  our  Equation  (1), 
but  with  the  cp  set  to  0,  which  would  tend  to  underestimate  the 
angular  extent  of  the  Ha-emitting  region.  The  other  possibility 
is  that  the  Ha-emitting  region  actually  grew  in  size  from  1994 
to  2007.  If  we  assume  that  the  difference  between  our  values 
and  those  reported  by  Vakili  et  al.  (1997)  to  be  representative 
of  a  real  change,  then  this  results  in  an  expansion  rate  of  the 
Ha-emitting  region  (due  to  changes  in  opacity)  over  this  13 
year  period  of  ~3  km  s-1.  This  is  well  within  the  expected 
wind  velocities  of  stars  with  radiation-driven  wind  structures, 
especially  P  Cyg,  which  has  a  terminal  wind  velocity  almost 
2  orders  of  magnitude  larger  (Barlow  et  al.  1994;  Lamers  et  al. 
1985).  Furthermore,  if  the  Ha-emitting  region  did  indeed  grow 
in  size,  we  would  expect  the  Ha  emission  to  be  stronger,  which 
is  supported  by  the  fact  that  the  Fmax/Fc  values  in  Table  3  that 
range  between  17  and  24  are  larger  than  the  ~  15  peak  value  of 
the  spectrum  taken  by  Vakili  et  al.  (1997). 
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4.2.  Signature  of  Asymmetry 

To  test  for  the  presence  of  a  signature  of  deviation  from  point 
symmetry  in  the  signal  in  the  Ha  channel,  we  have  plotted  the 
weighted  mean  closure  phases  in  Figure  3.  If  the  Ha  emission 
originates  from  an  intensity  distribution  that  is  not  symmetric 
across  the  origin  of  the  photocenter,  or  equivalently  the  Ha- 
emitting  region  is  not  concentric  with  the  central  star,  we  would 
expect  to  see  a  non-zero  closure  phase  in  the  Ha  channel.  On  a 
close  inspection  of  Figure  3,  we  conclude  that  the  Ha  closure 
phases  are  generally  very  close  to  0°.  Even  in  the  cases  where  the 
mean  closure  phase  deviates  from  zero  values,  it  is  only  at  the 
level  of  a  couple  of  degrees,  which  also  happens  to  be  only  two 
to  three  times  the  uncertainty  of  the  mean  value.  Therefore,  we 
conclude  that  based  on  our  observations  we  do  not  have  a  strong 
signature  of  deviation  from  point  symmetry  of  the  Ha -emitting 
region. 

Our  conclusions  are  not  necessarily  inconsistent  with  the  very 
similar  detection  of  a  subtle  phase  variation  (at  the  level  of  ~30°) 
across  the  Ha  emission  line  detected  by  Vakili  et  al.  (1997) 
who  interpreted  that  as  being  produced  by  a  localized  spatial 
feature  (i.e.,  a  localized  blob)  within  the  structure  of  the  wind. 
Assuming  that  the  closure  phase  variations  we  detect  across  the 
Ha  spectral  channel  in  Figure  3  are  caused  by  similar  spatial 
feature  (which  is  not  necessarily  expected  since  the  observations 
were  acquired  at  different  epochs),  we  would  indeed  expect  to 
see  a  much  weaker  signal  across  the  Ha  channel  for  two  reasons. 
The  first  being  that  we  measure  the  sums  of  three  phases,  which 
already  have  the  tendency  to  lower  the  signal  when  negative 
and  positive  phases  in  a  triangle  are  added  together.  The  second 
effect  is  related  to  our  much  wider  spectral  channel  that  contains 
the  entire  Ha  emission  line  and  a  significant  contribution  from 
the  central  star,  which  most  likely  can  be  described  well  by 
a  point-symmetric  intensity  distribution  that  contains  only  real 
components  in  the  Fourier  space  (i.e.,  only  real  phases).  The 
net  result  is  that  the  real  components  will  tend  to  lower  the  net 
complex  phase  detected  in  the  Ha  channel. 

Signatures  of  clumping  in  the  circumstellar  region  of  P  Cyg 
on  much  larger  spatial  scales  have  also  been  reported.  For  ex¬ 
ample,  based  on  near-diffraction-limited  observations  using  the 
AO  system  on  a  1.52  m  telescope,  Chesneau  et  al.  (2000)  re¬ 
ported  clumping  at  the  scales  of  200-600  mas.  Similarly,  P  Cyg 
was  also  imaged  using  the  Multielement  Radio-linked  Inter¬ 
ferometer  Network  (MERLIN)  array  by  Skinner  et  al.  (1997). 
Their  6  cm  observations  produced  images  of  the  circumstel¬ 
lar  structure  on  scales  of  100-200  mas  with  a  50  mas  resolu¬ 
tion,  revealing  structural  changes  on  a  40  day  timescale  along 
with  flux  variations  at  the  20%  level.  Skinner  et  al.  (1997)  ar¬ 
gued  that  the  structural  variations  could  not  be  attributed  to 
variations  in  the  mass-loss  rate.  They  suggested  recombination 
within  the  free-free-emitting  region  responsible  for  the  radio 
emission  could  explain  their  observations.  The  weak  signature 
of  the  non-zero  closure  phase  in  the  NPOI  data,  if  confirmed, 
could  indicate  that  the  dumpiness  observed  by  Chesneau  et  al. 
(2000)  and  Skinner  et  al.  (1997)  already  originates  on  scales  of 
less  than  10  mas  (less  than  25  stellar  radii). 

4.3.  Photometric  Variability 

Lastly,  to  assess  the  level  of  photometric  variability  during 
the  time  frame  covered  by  our  interferometric  observations, 
we  used  photoelectric  V-band  observations  from  the  American 
Association  of  Variable  Star  Observers.  The  V-band  light  curve 
from  2005  to  2009  is  shown  in  Figure  5  where  we  only  include 


data  for  which  the  measured  check  star  (HD  193369)  magnitude 
was  no  more  than  0.04  mag  away  from  the  mean  value  of 
5.573  mag.  We  see  both  a  mean  seasonal-timescale  variability 
on  the  order  of  ~0.02  mag  and  a  very  short-term  variability 
(for  observations  taken  within  days  to  weeks)  at  the  level  of 
±0.05  mag.  Although  the  short-term  variability  appears  to  have 
large  scatter,  the  photometric  quality  of  the  AAVSO  PEP  data  is 
generally  good,  with  internal  errors  on  the  order  of  5-10  mmag. 
Therefore,  based  on  the  photometric  data  we  conclude  that 
although  the  mean  photometric  level  appears  to  change  by  less 
than  0.02  mag,  short-term  variations  upward  of  0.1  mag  cannot 
be  ruled  out.  This  implies  that  the  continuum  level  in  P  Cyg  could 
vary  up  to  10%,  which  would  affect  directly  the  EW  measures 
at  the  same  level.  This  also  strongly  suggests  that  the  variations 
in  EW  measures  seen  in  Table  3  are  caused  by  variations  in  the 
continuum  level  and  not  the  emission  component  itself. 

5.  SUMMARY 

In  general,  mass-loss  rates  and  stellar  wind  properties  are 
extremely  important  characteristics  of  LB  Vs.  Mass-loss  rates 
for  P  Cyg  have  already  been  estimated  through  several  methods, 
but  none  of  themrely  on  the  spatially  resolved  Ha -emitting  inner 
region.  In  this  observational  study,  we  have  presented  the  most 
comprehensive  interferometric  data  set  on  P  Cyg  to  date.  Our 
results  provide  a  new  and  independent  source  of  observational 
constraints  that  can  be  utilized  when  determining  the  properties 
of  the  wind  structure  in  the  context  of  the  predicted  Ha  intensity 
distribution  on  the  sky. 

We  have  fitted  circularly  symmetric  UD,  Gaussian,  and  two- 
component  models  to  the  interferometric  observations  and  have 
shown  how  the  Ha-emitting  wind  structure  of  P  Cyg  cannot  be 
represented  by  a  simple  one-component  model.  We  conclude 
that  the  wind  structure  can  be  represented  fully  by  a  two- 
component  Gaussian  model  with  angular  diameters  of  0\  = 
5.64  ±  0.21  mas  and  02  —  1.80  ±  0.13  mas.  The  interferometric 
signature  can  also  be  equally  well  represented  with  a  two- 
component  model  consisting  of  a  UD  and  a  Gaussian  intensity 
distribution.  In  that  case,  the  angular  diameters  are  $UD  — 
3.06  ±  0.15  mas  and  0gd  =  5.46  ±  0.16  mas.  While  we 
cannot  conclude  if  a  double-Gaussian  model  or  a  UD  plus  a 
Gaussian  model  better  represents  the  wind  structure  of  P  Cyg,  it 
is  clear  that  both  models  represent  the  wind  structure  better  than 
a  single  Gaussian  or  a  UD  model.  We  conclude  that  P  Cyg’s 
wind  structure  is  complex,  possibly  containing  multiple  layers 
of  varying  opacities.  The  data  presented  in  this  study  might 
be  used  to  directly  constrain  wind  models  of  P  Cyg,  perhaps 
allowing  for  an  independent  observational  constraint  on  the 
/8-law. 

Based  on  our  spectroscopic  and  interferometric  observations, 
we  conclude  that  P  Cyg  is  stable  at  the  10%  level.  Based  on  the 
photometric  observations  we  also  rule  out  the  possibility  that 
large  changes  in  the  Ha  EW  measure  were  masked  by  changes  in 
the  continuum  level.  The  lack  of  variability  at  a  significant  level 
combined  with  the  spatially  resolved  Ha-emitting  region  that 
appears  to  be  stable  between  the  observing  seasons  is  suggestive 
of  an  extra-quiescent  phase  of  P  Cyg  between  2005  and  2008. 

The  Navy  Prototype  Optical  Interferometer  is  a  joint  project 
of  the  Naval  Research  Laboratory  and  the  U.S.  Naval  Obser¬ 
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by  the  Office  of  Naval  Research  and  the  Oceanographer  of  the 
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